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Abstract 

Saturation properties and liquid-gas phase transition of nucleus are analysed 
in the framework of Hatree-Fock theory. We modify Hill-Wheller formula with 
a finite-size-effect parameter by fitting the zero-temperature properties of nu- 
cleus. Employing Gogny effective interaction and phenomenological expression 
of Coulomb energy, we give the critical temperature of liquid-gas phase transi- 
tion of nucleus being about 12 MeV, which agrees with the result extracted from 
heavy- ion collision experiments. It is pointed out that a phenomenological for- 
mula of surface energy of hot nucleus is not available in the region where nucleon 
density is far away from the normal density. 
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I. Introduction 

The developments of heavy-ion pliysics have evoked great interest on the studies 
about statistic properties of infinite nuclear matter and nucleus, especially about equa- 
tion of state (EOS) and liquid-gas phase transition. There have been a lot of calcula- 
tions of infinite nuclear matter in the framework of Hatree-Fock theory (HFT), in which 
finite-size-effect and Coulomb interaction were ignored [1-10]. The finite-size-effect of 
nucleus and Coulomb interaction in nucleons should play important roles in realistic 



heavy-ion collisions. Refs. [0, |T2| took into account the finite-size-effect with the help 
of Hill-Wheller formula and Coulomb effect in the use of a phenomenological expression 
in the calculation of finite nuclear matter. It was found that the finite-size-effect leads 
to a reduction of critical temperature of liquid-gas transition in finite nuclear matter 
by about 5 ~ 10 MeV, and Coulomb effect leads to a reduction about 1 ~ 3 MeV. In 
our opinion, any approach which is applied to study hot nucleus though infinite nuclear 
matter approximation should be able to give reasonable zero-temperature saturation 
properties of nucleus also. Basing on the calculations of six typical nuclei (from 2oCa 
to ll^U) with several Skyrme effective interactions, Wang and Yang found that it 



is impossible to give correct bound energies and other saturation properties of nucleus 



by adopting Hill-Wheller formula directly as that in Refs. |11, Q2|. There is a general 



agreement that Gogny effective interaction is able to describe the long-range and 
medium-range behaviours of nucleon-nucleon interaction more reasonable than Skyrme 
interaction. In this paper, we will modified Hill-Wheller formula with a finite-size-effect 
parameter through the studies on the zero-temperature saturation properties of nucleus 
with Gogny interaction. Employing the modified Hill-Wheller formula, we study the 
liquid-gas transition of hot nucleus. We also analyse the applicability of a widely used 
phenomenological expression of surface energy of hot nucleus. 

The paper is organized as follows. In section II, we modified Hill-Wheller formula 
by fitting the zero-temperature properties of nucleus. We study the stability character 
and liquid-gas phase transition of hot nucleus in section HI. The surface energy of 
nucleus is discussed in section IV. As usual, the last section is reserved for summary. 

II. Formalism 

Let us begin with a brief review of Hatree-Fock theory of infinite nuclear matter 



system with Gogny interaction. Gogny Dl effective interaction is expressed as |jT5 

1=1,2 

+to (1 + XoP.) p"5(r) + tWLs{(Ti + fxs) ■ V X 5(r) y, (1) 

where p is nucleon density and P„, P^ are exchange operators of spin and isospin re- 
spectively. In Eq. (|l|), the part with 5 function is zero-range interaction which is 
density-dependent and the part with Gauss function is finite-range interaction which 
is density-independent. The parameters of Gogny Dl interaction are listed in table 1. 
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The Hatree-Fock single particle spectrum with Gogny Dl interaction is 



Ud 



+ Uq + Ud, 

2m 

J2{qq'\V\qq' -q'q)fq', 



dV 




dp 


qiqj - qjq^^ 



1 + exp [[3{tq + Uq + UD - p)] ' 



(2) 



where V is the volume of nuclear matter system, /i the chemical potential, T the 
temperature, /3 = 1/(A;bT), Uq the usual HF single particle potential, and ud the rear- 
rangement potential. It has been pointed out that the introduction of rearrangement 
potential is necessary to guarantee the density relation p = ^ J2q fq and give correct 
chemical potential Then EOS can be readily obtained from Eq. (j^) P, [12], |13|. 



Now we turn to the approach to respect finite-size-effect and Coulomb interaction in 
finite nuclear matter system. In Hatree-Fock theory, the calculation of infinite nuclear 
matter is simplified by taking plane wave function as single particle wave function, in 
which the state number between k and k + dk in the momentum space reads 



dN, = V 



edk 



(3) 



To study the liquid-gas phase transition in finite nuclear matter system, Jaqaman, 
Mekjian and Zamick [|TT| proposed that the finite-size-effect can be taken into account 
in the theory framework of infinite nuclear matter system by Hill-Wheller formula 



dNu. = V 



k'^dk _Skdk L dk 
~ V~8^ ^ V8^ 



(4) 



where S and L are the measures of the average surface and linear of finite nuclear matter 
system respectively. For a spherical system with radius R, S = AirR'^, L = 2ttR. The 
effects of zero-point motion and quantization of wave number in a finite-size system 
were taken into account in Eq. (^. In fact, Eq. is derived from a simple model, 
a finite-size system with free particles. Strictly speaking, it cannot be applied to the 
study of nucleus directly. The reason is as following: Because there are interactions in 
nucleons, the single particle wave function in nucleus is not plane wave function as that 
in the finite-size system without interaction, and different for the nucleus with different 
shell. Eq. (Q) is a approximate approach to study finite-size-effect in the theory 
framework of nuclear matter system. In our opinion, any approach which is applied 
to study the liquid-gas phase transition of nucleus through the infinite nuclear matter 
approximation should be able to give reasonable zero-temperature saturation properties 
of nucleus also. The calculations of the six typical nuclei show that Eq. can not 
satisfy this constrain (see table 2). It can be found that the bound energies per nucleon 
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of the nuclei calculated with Eq. (^) are much smaller than the experimental data and 
this divergence becomes more and more serious as the nucleon number decrease, so that 
2oCa can not combine any more. This shows that the finite-size-effect isn't considered 
suitably with Eq. (^). As a result, it is hardly to believe the calculations of EOS and 
liquid-gas transition of finite nuclear matter system given in Refs. [T^. We may 
add a finite-size-effect parameter ap to the second and third terms in Eq. (D to respect 
the fact that the single particle wave function in nucleus is different from that in the 
finite system with free particle. 



dm = V 



k^dk 



+ ap 



ATTp\ 3 3kdk 



3A 



in 



3 3dk 



IGtt 



(5) 



where ap is decided by fitting the experimental data of bound state energies of nuclei. 

Eq. @ becomes Eq. (|) when we take ap = I and becomes Eq. (|^) as A — > oo. 

By fitting the experimental data, we find that to take the finite-size-effect parameter 

ap = 0.35 in Eq. (|^), the bound energies of the nuclei from light to heavy agree with 

the experimental data very well and the saturation densities po of the six typical nuclei 

are about 0.13, which has mild dependence on the nucleon number A (see table 2). 

This agrees with the saturation character. With this po, through the phenomenological 

1 1 
formula R = r^A^ of the effective radius R = (3/47rpo)^, we obtain ro = 1.21 ~ 1.22, 

which agrees with the experimental result also. All of these show that it is reasonable 
with Eq. to study the finite-size-effect of nucleus. 

Coulomb interaction is a long-range interaction, adopting it directly in the cal- 
culation of matrix between plane wave functions will bring divergent result. So, for 
simplicity, we use a phenomenological expression of Coulomb energy per proton to 



respect Coulomb effect |11, 12 



c 



a Z' 



1-5 



WttZ 



pi A 



(6) 



We take a = 1.50 in order to fit the experimental data more well, which is slight 
different with the usual value a = 1.39. 



III. Stability character and liquid-gas phase transition of hot nucleus 

For discussing the stability character of hot nucleus, limit temperature T; is defined, 
which is the highest temperature below that free energy has a minimum. Nucleon 
density corresponding with the minimum of free energy at Ti is called limit density pi. 
When the temperature is larger than T;, the pressure of the system is always larger 
than zero. The condensation phase and the nucleon-gas phase can not arrive phase 
equilibrium without extra pressure. The values of Ti and pi of the six typical nuclei are 
listed in table 3. It can be found that Ti and pi are not sensitive to nucleon number A. 
Ti is about 9.0 ~ 9.5 MeV, and pi is about 0.075 fm~^. 

Liquid-gas phase transition of hot nucleus is another interesting question. Strictly 
speaking, a phase transition can only occur in a system with infinite number of particles. 
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which is reflected in the singularity behavior of some thermodynamic quantities. For 
example, the specific heat displays a sharp A-type singularity at the critical temperature 
Tc(oo) for a liquid-gas phase transition in an infinite particles system. The specific 
heat of finite particles system does not exhibit such a sharp singularity, but it has a 
large peak at a temperature Tc{A) which approaches to Tc{oo) as the particle number 
y4 — + oo. The temperature Tc{A) can be regard as the critical temperature of a finite 
particle system. For a finite particles system, it is convenient to determine the critical 
point by the inflection point condition of /i ~ p isothem rather than P ~ p isothem 



dp 



where p is the average chemical potential 



= 0, (7) 



1 1 — 7 1 + 7 

fJ'= -{PpfJ'P + PnPn) = ^—Pp + ^—Pn (8) 

with p denotes parton and n neutron. To study the finite-size-effect, we calculate 
the EOS and liquid-gas transition critical point of symmetrical infinite nuclear matter 
system with different nucleon number and without Coulomb interaction. Figs. 1 and 2 
display p ~ p isotherms at temperature T = 6 MeV and 14 MeV. As it is known, liquid 
phase and gas phase can coexist at low temperature, while only gas phase can exist at 
high temperature. It can be found that there is a critical nucleon number Ac{T) for a 
certain temperature. As A < Ac{T), only gas phase can exist. Table 4 lists the critical 
temperatures and densities pc of the symmetrical nuclear matter with different 
nucleon number A. and pc decrease with A decreasing. Here again we find that 
for the realistic nucleus system [A < 250) the liquid-gas critical points calculated with 
Eq. (^) and ap = 0.35 are very different from that with Eq. obviously. Coulomb 
interaction and neutron-proton asymmetry are another two factors which decrease the 
critical temperature. Table 3 gives the critical temperatures Tc and densities pc of 
the six typical nuclei. It can be found that with Eq. (^) and ap = 0.35, finite-size- 
effect. Coulomb interaction and asymmetry effect reduce the critical temperature 
by about 3 ~ 4 MeV comparing with that of the symmetrical infinite nuclear matter. 
Tc and pc are insensitive to nucleon number A, which is very different from the results 
with Eq. (p. As exhibited in section II, we can obtain reasonable zero-temperature 
saturation properties of nuclei with Eq. and ap = 0.35 while we cannot with Eq. 
(^. Thus the critical temperatures Tc obtained in the present work are reasonable and 
reliable. Tc has mild dependence on nucleon number A. In fact, all of the finite-size- 
effect. Coulomb interaction and neutron-proton asymmetry play the roles to weaken 
the nucleus combination and reduce Tc. The finite-size effect becomes weak with A 
increasing, while Coulomb interaction and neutron-proton asymmetry effect play more 
and more important roles. As a result, Tc is not sensitive to the nucleon number A, 
Tc ~ 12 MeV. Our results are close to that with Skyrme effective interaction WM. 
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Panagiotou et. al. |]T^ extracted the critical temperatures of nuclei by analyzing 



the mass distribution of multi-fragments in intermediate-energy heavy-ion collisions 
basing on the condensation theory, which gave T^. ~ 12 MeV. Bondorf et.al. []TB| 



obtained ~ 11 MeV for A=100 system basing on the Monte Carlo calculation. Our 
calculations are consistent with these results. 

IV. Surface energy of nucleus 

A phenomenological surface energy expression was applied to study the stable prop- 
erties and dynamical properties of nucleus in some works ]19|, pO[] , 

i;, = _ = 4.(-j -, (9) 

where a is the tension coefficient a ~ 1.2 MeV-fm~^ which is derived from the surface 
energy term in the bound energy phenomenological formula of nucleus. This approach 
takes the tension coefficient being independent of density, which gives the relation 

2 

Es ~ p~3. In the present calculations, the surface energy coming from the finite-size- 
effect can be derived naturally from the bound energy of nucleus and that of infinite 
symmetric nuclear matter, 

AE = E{A, p) - E{A ^ oo, p). (10) 



The surface energy of ^Zr obtained with Eqs. (y) and (|10|) are plotted in Fig. 3. 



Although Eqs. (^ and (IC) present similar results in the neighborhood of normal 
nucleon density Pq, their behaviors are very different in the other density regions. The 
saturation curves {E /A) ~ p, pressure-density curves P ^ p calculated with Eq. (|^) 
are plotted in Figs. 4 and 5, which are very different form that obtained form HFT. 
The saturation properties of nucleus calculated with Eq. (^ are given table 2. We find 
that according to Eq. (^, the bound energies agree with the experiment data, but the 
saturation densities of the nuclei are general larger than that of infinite nuclear matter 
and reduce as mass number A increasing, which is not reasonable obviously. All of 
these show that Eq. (P) is not available for the region where the nucleon density is far 
away from the normal nucleon density. We find that Eq. (|T0|) can be expressed by the 
following formula 



E's = Cipi A-^i + C2pi A-^ , (11) 

where Ci = 53.0 MeV-fm^ and C2 = 108 MeV-fm^ In Eq. (|ll]) the first term is 
derived from the fermion gas model and the finite-size-effect is taken into account by 
the second term. The density dependence of Eq. (pJ]) which reads E'^ ~ pa is very 
different form that of Eq. (|^). Tables 2 and 4 give the saturation properties and critical 
temperature in the use of Eq. (|Tl|), which agree with the results with Eq. (|]) very 
well. The calculations also show that the surface energy Eq. (|10D is not sensitive to 
temperature. Thus, Eq. (|1^) can be employed in the study about EOS of nucleus. 
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V. Summary and discussion 

Wc analyse the saturation properties and liquid-gas phase transition of nucleus and 
infinite nuclear matter in the framework of Hatree-Fock theory. At first, in order to 
give reasonable zero-temperature properties of nucleus, we modify Hill-Wheller formula 
which is used to respect the finite size effect of nucleus. Employing Gogny effective 
interaction and phenomenological Coulomb energy, we obtain the critical temperatures 
of liquid-gas phase transition of nucleus being about 12 MeV, which is consistent with 
the result extracted from heavy-ion collisions. The critical temperature of liquid-gas 
phase transition of nucleus has mild dependence on nucleon number A due to the 
cooperation of finite-size-effect. Coulomb interaction and neutron-proton asymmetry. 
In addition, the surface energy of nucleus is analysed. It is pointed out that a widely 
used phenomenological expression of surface energy is not available in studying hot 
nucleus. 
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i 


/ii(fm) 


Wi{MeV) 


Bi{MeV) 




Mi(MeV) 


1 


0.7 


-402.4 


-100 


-496.2 


-23.56 


2 


1.2 


-21.3 


-11.77 


37.24 


-68.81 


to = 1350 MeV-fm^ a = 1/3, Wis = 115 MeV-fm^, Xq = 1. 



Table 1. Parameters of Gogny Dl effective interaction. 









90 7». 




i«P6 


238 rr 
92 ^ 




Exp. 


8.55 


8.64 


8.71 


8.25 


7.87 


7.57 




Eq. d) 


8.97 


8.82 


8.72 


8.02 


7.45 


7.15 


Bound energy (MeV) 


Eq. ® 


-1.01 


-0.36 


0.66 


1.20 


1.30 


1.29 




Es 


8.77 


8.70 


8.70 


8.00 


7.41 


7.28 




E's 


8.76 


8.72 


8.73 


8.11 


7.56 


7.27 




Eq. d) 


0.131 


0.132 


0.134 


0.134 


0.134 


0.133 


Po(fm-3) 


Es 


0.185 


0.180 


0.177 


0.173 


0.170 


0.176 




E's 


0.129 


0.130 


0.133 


0.133 


0.133 


0.133 



Table 2. Zero-temperature saturation properties of the six typical finite nuclei. 



9 





40 


|M 


90 7», 






238 rr 
92 ^ 


SINM 


(1) 


T,{MeV) 


Eq. (1 
Eq. (I 




12.95 
5.10 


13.15 
6.35 


13.30 
7.55 


13.35 
8.50 


13.40 
9.05 


13.35 
9.20 


15.85 
15.85 


(1) 


Pc(fm-^) 


Eq. (1 
Eq. (1 




0.050 
0.025 


0.050 
0.030 


0.050 
0.035 


0.055 
0.040 


0.055 
0.040 


0.060 
0.040 


0.060 
0.060 


(2) 


T,{MeV) 


Eq. (1 
Eq. (1 




12.35 
4.35 


12.40 
5.35 


12.45 
6.50 


12.20 
7.30 


12.05 
7.60 


11.95 
7.65 




(2) 




Eq. (1 
Eq. (1 




0.050 
0.025 


0.050 
0.030 


0.050 
0.030 


0.050 
0.035 


0.050 
0.040 


0.050 
0.040 




(2) 


Ti{MeV) 
p;(fm-3) 


Eq. (1 
Eq. (1 




9.45 
0.075 


9.45 
0.075 


9.45 
0.075 


9.25 
0.075 


9.10 
0.075 


9.00 
0.075 





Table 3. The limit temperature T/, density pi and liquid-gas phase transition critical 
temperature T^, density pc of six typical nuclei. (1) without Coulomb interaction; 
(2) with Coulomb interaction. SINM=symmetric infinite nuclear matter. 



A 


50 


100 


200 


1000 


10000 


oo 




Eq. (D 


13.10 


13.55 


13.95 


14.65 


15.30 




%{MeV) 


E's 


12.65 


13.25 


13.80 


14.65 


15.30 


15.85 




Eq. (i 


5.15 


8.05 


9.70 


12.25 


14.20 






Eq. (i 


0.050 


0.055 


0.055 


0.055 


0.060 






E's 


0.050 


0.050 


0.055 


0.055 


0.060 


0.060 




Eq. (i 


0.025 


0.035 


0.040 


0.050 


0.055 





Table 4. Size dependence of critical point of liquid-gas transition in symmetric infinite 
nuclear matter. 
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Figure Captions 

Fig. 1 The /x ~ p isotherms of infinite symmetric nuclear matter at T = 6.0 MeV 
evaluated with nucleon number A = 100 (solid line), A = 200 (dashed line), 
A = 1000 (dotted line) and A = oo (dash-dotted line) respectively. 

Fig. 2 Similar as Fig. 1 and T = 14.0 MeV. 

Fig. 3 Surface energy of l^Zr at temperature T = 0.0 MeV evaluated with Eqs. ( [T0| ) 
and (D with ap = 0.35 (solid line); Eqs. and (D (dashed line); Eq. (|) 

(dotted line); and Eq. ( pH]) (dash-dotted line) respectively. 

Fig. 4 The saturation curves (E/A) ~ p calculated with Eq. (solid line), Eq. (^) 
(dashed fine), Eq. (^) (dotted line) and Eq. (|11]) (dash-dotted line) respectively. 

Fig. 5 The P ~ p isotherms. The explanation of the curves is similar to Fig. 4. 
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Fig. 2 
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Fig. 3 




Fig. 4 




Fig. 5 



